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Abstract 

We present the next-to-next-to-leading order post-Newtonian (PN) spin-orbit Hamiltonian for two 
self-gravitating spinning compact objects. If at least one of the objects is rapidly rotating, then the 
corresponding interaction is comparable in strength to a 3.5PN effect. The result in the present paper in 
fact completes the knowledge of the post-Newtonian Hamiltonian for binary spinning black holes up to 
and including 3.5PN. The Hamiltonian is checked via known results for the test-spin case and via the 
global Poincare algebra with the center-of-mass vector uniquely determined by an ansatz. 

PACS numbers: 04.25.Nx, 04.20.Fy, 04.25.-g, 97.80.-d, 45.50.Jf 

Keywords: post-Newtonian approximation; canonical formalism; approximation methods; equations of 
motion; binary stars 

1 Introduction 

In the present paper the next-to-next-to-leading order post-Newtonian (PN) spin-orbit Hamiltonian for two 
self-gravitating spinning compact objects is derived. This Hamiltonian is of the order 3.5PN if at least one 
of the objects is rapidly rotating. Indeed, the result in the present paper completes the knowledge of the 
post-Newtonian Hamiltonian (and thus of the equations of motion) for binary spinning black holes up to and 
including 3.5PN. Besides the well-known Newtonian and 1PN Hamiltonians, previous results to this order 
for point-masses are the 2PN [1-3], 2.5PN [4, 5], 3PN [6-10], and 3.5PN [11, 12] Hamiltonians. For the 
spin part the leading order can be found in [13-16] and the next-to-leading order in [17-19]. Atthe3.5PN 
level one also needs all Hamiltonians cubic in the spins derived in [20, 21]. Notice that so far these cubic 
Hamiltonians are known for binary black holes only, whereas all other mentioned Hamiltonians (including 
the one derived in the present paper) are valid for general compact objects [or have been generalized to this 
case, see [22-25] for the spin(l)-spin(l) level]. Further, tidal effects become very important for general 
compact objects like neutron stars, see, e.g., [26] and also [27]. 

The calculation in the present paper was performed within the canonical formalism of Arnowitt, Deser, 
and Misner (ADM) [28, 29], which was recently generalized to rotating objects at linear order in spin [30], 
see also [31-33]. For various other (noncanonical) derivations of post-Newtonian results at the point-mass 
level see [34-40] and references therein. The next-to-leading order in spin was also treated in [41^-7] . Also 
more than two compact objects have been treated at high post-Newtonian orders [48-52]. Notice that the 
calculation in the present paper is comparable in complexity to the one of the 3PN point-mass Hamiltonian. 
In particular one has to check for certain integral contributions that can only be handled correctly within 
dimensional regularization, as they may lead to ambiguities when integrated in three spatial dimensions 
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[10, 53]. The result of the present paper was verified using the test-spin Hamiltonian in Kerr-spacetime 
given in [54] and using the global Poincare algebra. For the latter the center-of-mass vector had to be 
determined uniquely from an ansatz with 68 coefficients. 

It should be noted that the Hamiltonian derived in the present paper is yet only useful within the Taylor- 
expanded post-Newtonian series if at least one of the objects is rapidly rotating (due to the missing 4PN 
point-mass Hamiltonian). However, if the 4PN point-mass Hamiltonian can be derived in the future, then 
the result of the present paper must be included in the post-Newtonian series also when the spins are small. 
But the effective one-body approach for nonspinning objects is able to cover such higher post-Newtonian 
orders by calibration to numerical relativity, see [55, 56] and references therein. The result given in the 
present paper is thus expected to be valuable for the effective one-body formalism, just like the next-to- 
leading order one [57-59]. 

The paper is organized as follows. In Sect. 2 a brief outline of the calculation is given. The next-to- 
next-to-leading order spin-orbit Hamiltonian is presented in Sect. 3. A comparison with known results for 
the test-spin case is made. The Hamiltonian is further checked via the global Poincare algebra in Sect. 4, 
where the center-of-mass vector is uniquely determined from an ansatz. In forthcoming papers we will 
derive the next-to-next-to-leading order spin(l)-spin(2) Hamiltonian and provide much more details on the 
calculation of the spin-orbit Hamiltonian, shown in the present paper, as well. 

Three-dimensional vectors are written in boldface and their components are denoted by Latin indices. 
The scalar product between two vectors a and b is denoted by (ab) = (a • b). Our units are such that 
c = 1. There is no special convention for Newton's gravitational constant G. In the results P Q denotes the 
canonical linear momentum of the ath object, z a the position of the object, m a the mass of the object, S a 
the spin of the object, r a b = \z a — Zb\ the relative distance between two objects, and n a b = (z a — Zb)/r a b 
the direction vector pointing from object b to object a. In the integrands r a = |x — z a |, n a = (x — z a )/r a , 
and s a b — r a + rb + r a b- In the binary case the object labels a, b take only the values 1 and 2. 

2 Outline of the calculation 

In the following we present a short outline of our calculation - which will be discussed more detailed in a 
forthcoming publication - and cite the main literature necessary to undertake it. 

For all computations we used xTENSOR [60], a free package for MATHEMATICA [61], especially 
because of its fast index canonicalizer based on the package xPerm [62]. We also used the package XPERT 
[63], which is part of xTENSOR, for performing the perturbative part of our calculations. Furthermore we 
wrote several MATHEMATICA packages ourselves for evaluating integrals. 

First we generalized the derivation of the canonical formalism given in [30] to arbitrary dimensions. 
The initial action is of the same form as in d = 3 and the (d + 1) -split is also straightforward. The d- 
dependence enters the calculation via the relation between extrinsic curvature and field momentum, and 
via the decomposition of the metric and field momentum in the ADM transverse-traceless gauge. The 
calculation was done in d dimensions due to possible appearance of ambiguities in three-dimensional in- 
tegrals. Ambiguity means that one will get different results when one does an integration by parts in a 
certain integral. These ambiguities can only be ruled out or corrected by doing the UV-analysis explained 
in [10, 53], which relies on the complete d-dependence of the integrands. In the following UV-analysis 
always refers to the analysis of certain integrals via a Taylor expansion of the field expressions in the r\ 
variable, extracting the critical T\ powers, and averaging over the ni vectors afterwards to get the pole part 
of the integrals in e = d — 3. 

The integrations by parts necessary to get the Hamiltonian presented in this paper were done like 
suggested in [6, 10, 32] to get comparable intermediate results. After accomplishing the integration by 
parts we ended up with a Hamiltonian density which can be split up into three parts: A kinetic field part 
(containing the kinetic energy of the propagating field degrees of freedom) a matter part (containing only 
field-matter-interactions), and an interaction part (containing interactions between matter fields and the 
propagating fields), see, e.g., [4, 5]. We checked intermediate results against [6, 30, 32, 33]. 

From the Hamiltonian obtained by integrating the density from above one can go to a Routhian (a 
Hamiltonian in matter degrees of freedom and a Lagrangian in the propagating field degrees of freedom) as 
suggested in [6, 33] and can eliminate the propagating degrees of freedom by inserting their approximate 
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solutions in terms of the matter variables [6]. Subsequent elimination of time derivatives corresponds to a 
coordinate transformation [64, 65]. 

After obtaining a suitable Hamiltonian density via the simplifications mentioned in the last paragraphs, 
one has to integrate all appearing terms in the density to get a full Hamiltonian. The appearing integrals can 
be divided into three types: the delta-type J d d xf(x)Si, the Riesz-type J d^xn^ . . . n^n? 2 ■ ■ ■ n^rfr^ 
and the generalized Riesz-type J d 3 x n 1 ^ . . . n^n^ 1 . . . n J 2 l r" r 2 sj 2 . 

The delta-type integrals can be solved by the Partie-Finie regularization procedure mentioned in the 
appendices of, e.g., [6, 11, 33]. Another possibility to solve some of them is via the Riesz kernel method, 
where one inserts a Riesz kernel for the delta functions [53, 66]. This was done as an alternative way to 
check whether our algorithms work correctly, although we did not calculate all delta-type integrals via the 
Riesz kernel method because not all appearing inverse Laplacians can be solved via the method mentioned 
in [ 1 1 ] . Of course, using the Riesz kernel instead of a delta as source of the fields (to eliminate the necessity 
of distributional derivatives) makes the integration much more complicated because all delta-type integrals 
will be changed into integrals of the Riesz-type or the generalized Riesz-type. 

The Riesz-type integrals can be solved by eliminating the ni and n 2 vectors via rewriting them into 
derivatives as shown in [6, 33] and solving the remaining scalar integrals via the Riesz-formula [53, 66] 
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The ni and n 2 vectors in the integrands of generalized Riesz-type cannot be eliminated via rewriting the 
vectors into derivatives. Instead one has to use the averaging procedure in prolate spheroidal coordinates 
in [6] to get rid of the n vectors. Afterwards one can use the generalized Riesz formula which was found 
by R Jaranowski during his 3PN point-mass calculations (also given in [6]), 
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which reduces to the formula for the integrals of the Riesz-type for 7 — > 0. The function Ii/2(x, y) is the 
regularized incomplete Euler beta function which is defined as 

h/2(x,y)= R , , , (3) 



B{x,y) 



with 



B 1/2 (x,y) = ±. 2 F 1 (l-y,x,x + l;^, B(x, y) = , (4 ) 

being the incomplete Euler beta function and the Euler beta function respectively. Notice that all integrals 
of the generalized Riesz-type and all inverse Laplacians of two variables can only be solved in d = 3. Only 
the UV-singular part of those integrals can be evaluated in d dimensions. 

It turns out after using the integration procedures mentioned above that all integrands of the generalized 
Riesz-type appearing at spin-orbit level have such a structure that the incomplete Euler beta functions 
appearing there can be reduced to Gamma functions and Polygamma functions, which could be handled 
very well by Mathematica. 

3 Result 

To check our code we recalculated parts of the 3PN point-mass Hamiltonian given in [6, 10]. The next-to- 
next-to-leading order spin-orbit Hamiltonian we obtained as a result of the procedures discussed in Sect. 2 
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Obviously there are no logarithmic dependencies of ria appearing. 1 Also notice that the canonical anti- 
symmetric spin tensor was rewritten in terms of the canonical spin vector, which is possible in d — 3. 
The d-dimensional UV-analysis described in [10, 53] and Sect. 2 gave contributions to intermediate ex- 
pressions, however they exactly canceled in the final result. In contrast, for point-masses only the poles 
in e = d — 3 canceled but a finite part remained. Further from a combinatorial point of view there are 
66 algebraically different possible contributions to the Hamiltonian for each object (written in terms of 
the canonical spin tensor), but 24 of them do not appear in the canonical representation used here. The 
Hamiltonian is valid for any compact objects like black holes or neutron stars. It completes the knowledge 
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of the dynamics up to and including 3.5PN for maximally rotating black holes. For other objects the S 3 
Hamiltonians and the inclusion of tidal effects are missing. One can find a discussion of leading order tidal 
effects in [26]. Notice that the coupling structure in this Hamiltonian reduces in the center-of-mass frame 
to a pure (L S) structure with complicated coefficients. So the Hamiltonian is indeed a spin-orbit Hamilto- 
nian. We compared our result in the test-spin limit with the PN expanded Hamiltonians of a test-spin near 
a Kerr black hole in ADM coordinates in [54, Eq. (6.20)] and got full agreement. 2 

The matter variables appearing in the fully reduced matter- only Hamiltonian fulfill the standard Poisson 
bracket relations, namely 

{zliPaj} = , {S a (i), S a (j)} — SijkS a (k) , (6) 

all other zero and the Hamiltonian can be used to get the time evolution of an arbitrary phase space function 
A via 

Although the algorithms to be used at formal 3PN level are given above, the whole calculation is very 
hard. It is in particular much harder than (and really different from) the formal 2PN calculation at next-to- 
leading order spin-orbit level. 



4 Approximate Poincare algebra 



In this section we check that the global Poincare algebra is fulfilled in a PN approximate way, see, e.g., 
[9, 17]. Besides the Hamiltonian, the quantities entering the Poincare algebra are the center-of-mass vector 
G, the total linear momentum P, and the total angular momentum tensor J y = — J Ji . As the latter two 
are the infinitesimal generators of translations and rotations, they can be expressed in terms of canonical 
variables as [9, 17] 



P = ^P a , J* 



(8) 



where the canonical spin tensor S a u)(fi is related to the canonical spin vector S a via S a = EijkS a (m. 
Notice that the total angular momentum is not only the sum of the orbital angular momenta but also contains 
the spin angular momenta. For the contributions of the propagating field degrees of freedom see, e.g., [31, 
32]. As in [9, 17] we used an ansatz for the center-of-mass vector G at next-to-next-to-leading order spin- 
orbit level which contains 68 unknown coefficients here. For comparison we mention that the 2PN binary 
point-mass G-vector requires 20 unknown coefficients to be fixed and the 3PN binary point-mass G-vector 
requires 78 unknown coefficients to be fixed [9]. In [17, Eq. (5.9)] one can see the ansatz for the next-to- 
leading order case (which contains only 8 unknown coefficients). But at the order considered here there will 
be additional linear momentum powers which increase the number of necessary coefficients significantly. 
16 of them can be fixed by taking into account the {G*,P^} Poisson bracket relation appearing in the 
Poincare algebra. The remaining 52 coefficients were uniquely fixed by evaluating the {G, H} Poisson 
brackets. The consistency of the solution obtained by evaluating the Poisson bracket above was checked 
by evaluating the {G 1 , G J } Poisson bracket relation and all remaining relations of the Poincare algebra. 
The center-of-mass vector at next-to-next-to-leading order spin-orbit level is given by 
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